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ABSTRACT 

We present two weak lensing reconstructions of the nearby {zd = 0.055) merging cluster Abell 
3667, based on observations taken ^ 1 year apart under different seeing conditions. This is the 
lowest redshift cluster with a weak lensing mass reconstruction to date. The reproducibility of 
features in the two mass maps demonstrate that weak lensing studies of low-redshift clusters are 
feasible. These data constitute the first results from an X-ray luminosity-selected weak lensing 
survey of 19 low-redshift (z < 0.1) southern clusters. 

Subject headings: galaxies: clusters: individual (Abell 3667), gravitational lensing 



1. Introduction 

Weak gravitational lensing — measurement of 
the induced shear of distant galaxy images to 
infer the foreground mass distribution — has now 
become a standard tool to probe dark matter in 
galaxy clusters. Since the work of Tyson et al. 
(1990), weak lensing studies of over twenty galaxy 
clusters have been pubhshed (Melher 1999). In all 
cases to date, the observations have been confined 
to moderate to high-redshift clusters, Zd > 0.15*, 
for two reasons: (i) the cluster lensing strength 
is maximized if the angular diameter distance to 
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the foreground cluster is roughly half that of the 
background source galaxies, which at a limiting 
depth of i? ^ 25 are typically at (zg) ^ 1; (ii) a 
distant cluster subtends an area on the sky which 
can be roughly encompassed by a single CCD 
chip, allowing deep exposures of the entire field 
in moderate observing time. Unfortunately, most 
of our knowledge of cluster properties comes from 
X-ray and dynamical (optical redshift) studies, 
which are much more easily conducted for nearby 
clusters. This mismatch between the observing 
requirements for weak lensing and dynamical 
studies has hampered direct comparison of these 
mass measuring techniques. It would be beneficial 
to have weak lensing observations for a sample 
of nearby clusters. This would enable more 
detailed studies of the cluster binding mass, 
baryon fraction, and morphology. 

For weak lensing studies, nearby clusters have 
the advantage that their background galaxies are 
relatively well resolved (reducing the effect of PSF 
smearing) and their inferred projected mass distri- 
butions are relatively insensitive to the uncertain 
background source galaxy redshift distribution in 
the limit Zd ^ Zg. The advent of CCD mosaic 
cameras on 4m telescopes, coupled with develop- 
ments in weak lensing analysis, now make possible 
wide-field studies of weak lensing in nearby clus- 
ters (Joffre et al. 1999; Stebbins et al. 1996). 

With these benefits in mind, we have begun 
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a complete, X-ray luminosity-selected weak Icns- 
ing survey of 19 nearby {zd < 0.1) Southern 
clusters (Joffre et al. 1999); this sample, drawn 
from the XBACS (Ebehng et al. 1996), is also 
being targeted by the Viper telescope for Sunyaev- 
Zel'dovich observations (Romer et al. 1999). As 
of October 1999, we have completed imaging ob- 
servations of 9 clusters in the sample, using the 
BTC and Mosaic II cameras at the 4m telescope 
at CTIO. In this Letter, we present results for our 
first target, the z = 0.055 cluster Abell 3667. 

The major obstacle to weak lensing studies of 
nearby clusters is the small shear signal they pro- 
duce. Shear estimates are limited by the intrinsic 
ellipticity of the background galaxies (shape noise) 
and PSF anisotropics. To create reliable mass 
maps we must therefore measure ellipticities for a 
large number of background galaxies and correct 
for the PSF anisotropy vc;ry accurately. Given 
that we expect cluster-induced shears of at most ~ 
5%, we have reduced systematic sources of shear, 
which include atmospheric turbulence, telescope 
shake, and optical aberration, to less than 0.5% 
across the images (based on the corrected stellar 
ellipticities). As proof of our ability to make these 
measurements, we present weak lensing analyses 
of A3667 based on two sets of observations taken 
~ 1 year apart under different seeing conditions 
and with different observing strategies. 

2. Observations 

Much ancillary information is available for 
Abell 3667. Based on redshifts for 154 cluster 
members (Sodre et al. 1992; Katgert et al. 1996), 
an optical velocity dispersion in the range fT„ = 
970-1200 km/s has been inferred (Sodre et al. 1992; 
Fadda et al. 1996). From the ROSAT All-Sky 
Survey, the X-ray luminosity Lx = 6.5 x 10^^ 
ergs/s (Ebeling et al. 1996), among the strongest 
X-ray sources in the southern sky. The radio 
map shows a double halo structure (Rottgering 
et al. 1997). The A3667 field contains two 
dominant D galaxies; the ROSAT PSPC im- 
age (Knopp et al. 1996) and ASCA temperature 
map (Markevitch et al. 1998) indicate the clus- 
ter is undergoing a merger along the direction 
(Knopp et al. 1996) joining them. 

We observed A3667 on two separate nms with 
the BTC (Wittman et al. 1998) at the CTIO 4m 



telescope. The BTC is a 4096 x 4096 pixel camera 
with a pixel scale of 0.43"; there are significant 
gaps between the four chips that must be re- 
moved from the final combined image by dithering 
exposures. In June 1997, two of us (PF and 
JM) observed A3667 in the R and Bj bands in 
relatively poor seeing (combined stellar FWHM 
of 1.55" ± 0.15" in R). The combined image in 
each filter covers an approximate area 42' x 42', 
with a maximum surface brightness limit of 28.6 
mag arcsec"^ in R and 28.9 mag arcsec"^ in Bj, 
corresponding to Icr in the sky. Total observing 
times were 5500s in Bj and 15500s in R. We 
denote these the 'a' set of images in our analysis. 
In September 1998, three of us (JF, TM, RN) 
observed A3667 under better seeing conditions 
(combined stellar FWHM = 1.23" ± 0.07" in i?), 
covering a 44' x 44' area to a maximum surface 
brightness limit of 28.2 mag arcsec"^ in R and 
27.9 mag arcsec^^ in B. This corresponded to 
a total observing time of 12600s in R and 2250s 
in B. We also obtained /-band images which we 
have used to construct a color-magnitude diagram 
to remove cluster members from the background 
sample. (Due to fringing effects, the /-band data 
was not used in the lensing analysis itself.) We 
denote these better-seeing images the set. 

After debiasing and flatfielding of the frames, 
the frames were coadded. BTC observations 
of USNO astrometric fields were used to re- 
move field distortion from the images. Ob- 
jects were detected using SExtractor (v2.1.0) 
(Bertin & Arnouts 1996). Our own code measures 
quadrupole moments Qij of each object using the 
adaptive Gaussian weighting scheme of Bernstein, 
et al.(2000): the image is multiplied by an 
elliptical Gaussian weight function, the shape 
of which is successively iterated to match the 
object's ellipticity. This routine returns estimates 
of the ellipticity vectors ei = (Qn — (522)/(Qii + 
(^22) and 62 = 2(3i2/((3ii + Q22) and their 
uncertainties. The unsaturated bright stars in 
each imago were used to characterize the PSF 
anisotropy as a function of position. Following 
Fischer & Tyson (1997), we convolved the images 
with a spatially varying kernel which circularizes 
the stellar images. We then repeat the detection 
and measurement of background galaxies in the 
PSF-corrected images. The moments of each 
galaxy are finally corrected for isotropic PSF 
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dilution using the simulations and analytic results 
of Bernstein, et al.(2000); the correction factor 
depends upon galaxy image size relative to the 
PSF and upon the image profile. Objects in the 
background galaxy samples for the lensing analysis 
are selected to have magnitudes R < 24.75, 
B < 24.5, and Bj < 25.5. They arc also required 
to have half-light radii at least 1.5 times that of 
the PSF. The R — I color-magnitude diagram 
is used to remove red chistcr members from 
the background sample brighter than R = 22. 
These cuts ensure that the object moments can 
be accurately measured {S/N > 8 — lOcr), that 
stellar contamination is minimal, and that the vast 
majority of the sample lies well behind the cluster. 
The resulting samples contain approximately 
11,000 background galaxies for the a R image, 
30,000 for the /3 R, 18,000 for the a Bj, and 
11,000 objects for the P B image. We note that 
the (3 R image has a larger number of galaxies 
than the a set, although the later has a longer 
exposure time. This is due to the fact that the 
P images covered a larger area on the sky and 
the poor seeing of the a R image makes galaxies 
and stars difficult to distinguish. Many of these 
smaller galaxies are lost when we cut on half-light 
radius. 

3. Comparison 

To study the robustness of the mass maps 
derived from the two sets of observations, we 
first examine the consistency of the measured 
ellipticities of background galaxies. We focus on 
the shear measurements (as opposed to the mass 
maps themselves) because we expect the errors 
to be uncorrelated in different regions of the sky. 
In principle, differences in seeing, imaging depth, 
and filters will lead to differences between the 
shear fields estimated from the two observations. 
However, after applying the corrections discussed 
above, the derived shear fields should be strongly 
correlated. 

We trimmed the R and B band images in the a 
and P datasets to the regions common to both and 
calculated the mean ellipticities of the background 
galaxies in 100 angular bins of width ^ 25". 
To quantify the consistency of the two fields we 
calculate the value for the ellipticities for each 
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The estimate of the ellipticity variance in a spatial 
bin for a given data set is 

al, = {NjN^){al), + (Na/Nfal^J{Na - 1), 

where N is the total number of background galax- 
ies in the bin and Nc and Nd are the numbers 
of galaxies in the bin with and without coun- 
terparts in the other data set. (o'm)c is the 
average measurement error in derived from our 
measurement uncertainties in Qij for the galaxies 
common to both sets; fr^^s is the rms spread in 
the ellipticities of the galaxies in the bin not 
found in the other set. That is, for galaxies mea- 
sured in both data sets, we use the measurement 
uncertainty rather than the rms per bin to take 
into account the correlation between the data sets. 

Figure la shows the binned ellipticities of back- 
ground galaxies found in both data sets in the 
R band; the reduced is l-H for 200 d.o.f. 
Fig. lb compares the ellipticities for all galaxies: 
in this case, the scatter between the datasets is 
significantly larger, as expected since they are no 
longer confined to the same galaxies; the reduced 
X^ for all the objects is 1.19, implying a probability 
P(X^|200) = 3.5% in the case of random errors. 
The average ellipticity in Fig. lb is 2.57 ± 0.22% 
for the a set and 2.42 ± 0.16% for the /? set. 
As expected for low rcdshift cluster lensing, the 
maximum tangential ellipticity is only 4%, while 
the majority of the signal is ^ 1%. The major 
contribution to x^ comes from a few bins located 
at the edge of the images; removing the two worst 
bins drops substantially, raising the probability 
that the ellipticities are consistent within the er- 
rors to ~ 30%. The edges of our images have the 
largest field distortions and shallowest coverage 
and are therefore where we expect the largest 
discrepancy between the data sets. For the blue 
filters, due to the shallowness of the /3 set and the 
fact that the B and Bj filters are not identical, 
there are very few background objects in common 
(~ 1000 vs. 9000 in the R images). The error 
in x^ ii^ dominated by the rms ellipticity, giving 
a value of x^ = 1-125 or P{x^\2m) = 10.9%. 
This is a very strict test of our measurement 
and correction algorithm, showing that our shear 
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measurements are reprodueible. The relatively 
low formal probabilities may be due to either small 
residual systematics or slight underestimates of 
the measurement error. Across the image, the 
high correlation between the corrected data sets is 
quite encouraging, with large differences generally 
confined to small areas near the edge of the field. 

We performed mass reconstruction on a 60 x 
60 grid, applying a version of the Kaiser and 
Squires (Squires et al. 1996) algorithm separately 
to the blue and red catalogs and combining them 
by weighting the mass at each gridpoint by its 
S/N. We have chosen to plot the S/N of each 
mass map pixel as this gives a direct picture 
of which mass peaks are significant. We esti- 
mated the the noise in each pixel for this map 
as follows: we rotated each galaxy orientation 
through an arbitrary angle and then computed 
the resulting mass map; we then repeated this 
procedure 100 times and estimated the noise from 
the variance of these 100 noise maps. As a 
check of systematics, wc produced mass maps 
for each image with all galaxy orientations ro- 
tated by 45 degrees and found no significant fea- 
tures. In the absence of residual systematic ef- 
fects such a map should be consistent with noise 
(Stebbins et al. 1996; Kaiser et al. 1994). As an 
additional systematics check, we have made shear 
maps for the stars in each filter; they are consistent 
with noise, as expected. 

The combined convergence maps are shown in 
Figures 2 and 3 for the a and (3 data sets. There 
is remarkable agreement in morphology of the two 
maps, and the largest mass features appear to be 
robust. To quantify this correlation, wc calculate 
Pearson's r between the two mass maps: 

^ _ Y^^jnia - {ma)){mi3 - {mp)) 

where m is the value of a mass map pixel and 
N is the total number of gridpoints. r = 
— 1 for completely anticorrelated data, for 
uncorrelated data, and 1 for completely correlated 
data. We find a value of r = 0.6 between 
the two maps, indicating that they are in fact 
correlated. The formal probability of achieving 
such a high value of r by chance for uncorrelated 
maps is erfc(|r| y^N/2), negligibly small for 3600 
gridpoints. To further quantify the degree of 



correlation, we calculated Pearson's r between 
the 100 noise maps described above and found 
these maps gave a value of —0.009 ± 0.089. To 
determine the correlations between the entire 
maps, rather than the correlation introduced by 
the coincidence of the central peak, we masked 
out a box centered on the two maps with a 
size of 10'. When Pearson's r was computed on 
the remaining unmasked areas, we still found a 
value of r — 0.40. This value remained fairly 
constant as we increased the masked region until 
the masked region's size approached that of the 
two images. Inspection of the maps indicates 
that even the relatively low significance (~ 3cr) 
mass peaks which correspond to a convergence 
of ~ 0.02 are reproducible. We quantified this by 
masking regions greater than 3a in either map and 
computing Pearson's r on the unmasked portions. 
The value of the correlation only dropped to 
r = 0.27. We note that r dropped quickly to zero 
if we masked out all regions of significance lower 
than 3cr. We have also implemented the mass 
reconstruction algorithm of Seitz & Schneider 
(1998) and find the same correlation between 
features. 

The mass map is strongly peaked around the 
central D galaxy and generally correlates well with 
both the cluster light and ROSAT X-ray flux 
distributions (Joffre et al. 2000). At lower signif- 
icance, there also appears to be mass associated 
with the second bright D galaxy in the NW of the 
image and with cluster galaxies in the N and SE 
of the central D. 

4. Conclusion 

We have detected weak lensing at the 0.5 — 4% 
level in a nearby galaxy cluster using two sepa- 
rate sets of images. Despite differences in depth, 
seeing, and filters, the shear maps are consistent 
within the errors, and the reconstructed mass 
maps are strongly correlated. This reproducibility 
demonstrates the feasibility of using weak lensing 
to probe the mass distribution in low-redshift 
clusters. In future work, we will apply these 
methods to all the clusters in our survey and 
compare the resulting mass maps with X-ray, 
optical, and Sunyaev-Zel'dovich data to obtain a 
more detailed picture of the properties of nearby 
clusters. 
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Fig. 1. — (a) The binned average ellipticities of background galaxies appearing in both R band images, light 
blue for the a data set and black for the /3 data set. (b) The same for all background galaxies. The units 
are CCD pixels; an ellipticity of 2% is shown in the lower left. 
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Fig. 2. — Projected mass map of A3667 for set a. Contours correspond to a change in signal to noise of 
unity; only contours > la are shown. The mass map is superimposed upon the a R image. The image is 
42' X 42'. 
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Fig. 3. — Projected mass map of A3667 for set [3. The mass map is superimposed upon the [3 R image. The 
image is 44' x 44'. 
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